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^vq We study some problems related to the effect of bounded, additive sample noise in the ban- 

dlimited interpolation given by the Whittaker-Shannon-Kotelnikov (WSK) sampling formula. 
We establish a generalized form of the WSK series that allows us to consider the bandlimited 

■ \ interpolation of any bounded sequence at the zeros of a sine-type function. The main result 

of the paper is that if the samples in this series consist of independent, uniformly distributed 
random variables, then the resulting bandlimited interpolation almost surely has a bounded 
global average. In this context, we also explore the related notion of a bandlimited function 
with bounded mean oscillation. We prove some properties of such functions, and in particular, 
we show that they are either bounded or have unbounded samples at any positive sampling 
rate. We also discuss a few concrete examples of functions that demonstrate these properties. 
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1 Introduction 



The classical Whittaker-Shannon-Kotelnikov (WSK) sampling theorem is a central result in signal 
processing and forms the basis of analog-to-digital and digital-to-analog conversion in a variety 
of contexts involving signal encoding, transmission and detection. If we normalize the Fourier 
transform as /(ft)) = J^ 00 f(t)e~ 2nmt dt, then the sampling theorem states that a function / e L 2 (IR) 
with supp(/) C [— |, |] can be expressed as a series of the form 

,, , f, sm(n(bt-k)) 

m= k L ak «(»-*) ' () 

/2 



where a^ = f(k) are its samples. Conversely, for a given collection of data {a{\ G I , the series dl 

b b 

2' 2 



defines a function in L 2 (M) with supp(/) C [— |, |] called the bandlimited interpolation of {a^}. 
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The calculation or approximation of this series is a standard procedure in many applications. For 
example, in audio processing it is used for resampling signals at a higher rate, typically by applying 
a lowpass filter to the piecewise-constant zero order hold function of the samples [6]. In this paper, 
we consider the situation of bounded noise in the samples a&. Building on recent work by Boche 
and Monich on related problems E HI 13, we study some properties of the effect of the noise on 
the bandlimited interpolation /. 



Before we discuss our problems, it will be convenient to define the Paley-Wiener spaces for 1 < 
p < oo by 

PW*= {/ 6 LP : supp(/) C - h - b - 

where / is interpreted in the sense of tempered distributions. Our notation FWE essentially fol- 
lows Seip lTT2ll . and is slightly different from the one used by Boche and Monich. Without loss of 
generality, we will set b = 1 in what follows. 

Returning to the series ([[]), we consider corrupted samples of the form a\ = Ti +Nk, where 
Tk are the true samples and A/)t is some form of noise, and we correspondingly write f(t) = 
T(t) +N(t). One obstacle we face is that the noise {A/*} may not naturally decay in time alongside 
the signal, and even if {7^} G I 2 , it is often more physically meaningful to consider {A^} G l°°. The 
WSK sampling theorem shows that for any collection of samples {a^} G I 2 , there exists a unique 
function / G PW 2 with f[k) = a&. However, for bounded samples {a^\ G l°°, the series <TTj> does not 
necessarily converge. In fact, a given {a^} G l°° may correspond to multiple functions / G PW™, 
or to no such function Q. 

A simple example of the former possibility (non-uniqueness) is given by a^ = 0, which cor- 
responds to the functions fit) = and f[t) = s'm(nt). It turns out that adding one extra sample 
to the collection {a^} resolves this ambiguity, and allows us to consider the unique bandlimited 
interpolation of any bounded data {a^} G l°°. We discuss the details of this procedure in Section 3. 
The latter possibility (non-existence) is less obvious, but in [3], Boche and Monich presented an 
explicit example of this phenomenon. They showed that for the samples given by % = 0, k < 1, 
and at = ( — l)*yiog(fc4- 1), k > 1, there is no / G PW™ with f[k) = a^. It is also possible to con- 
struct other, similar examples using standard special functions, and we describe one such sequence 
of {a/c} in Section[3]and discuss its properties. 

The main observation of this paper is that such examples of {%} are in a sense "highly oscil- 
lating." By assuming that the noise N^ is statistically incoherent and defining N(t) carefully, we 
can rule out these examples and obtain sharper statements on the behavior of N(t). More precisely, 
we show in Section [4] that if Nk is a uniformly distributed, independent white noise process, then 
sup,. >0 ^ \ r _ r \N{t) \dt < °° almost surely. In other words, the average of \N(t) \ is globally bounded. 
We find that this result does not generally hold for {N^} G l°° that lack such a statistical condition, 
and we discuss examples that illustrate the differences. 

We also study a second topic motivated by further understanding N(t). As discussed in Q, 
the WSK series ([TJ) can be interpreted as a discrete Hilbert transform operator H, mapping a space 



of samples into a space of bandlimited functions (see also Q] and IfTTTO . The Plancherel formula 
shows that H maps I 2 into PW\ . In fact, H also maps l p into PW[ for any 1 < p < °«, and the series 
([T]) converges for any {a^} G / p [fTOll . This can be compared with the continuous Hilbert transform, 
and more generally any Calderon-Zygmund singular integral operator, which maps LP into itself 
for any 1 < p < °°. Such operators behave differently for p = °°, mapping L°° into the space BMO 
of functions with bounded mean oscillation lfT3ll . 

It is thus reasonable to expect that if we consider samples {a^} G l°°, the "right" target space for 
H may be one of bandlimited functions lying in the space BMO. However, this heuristic reasoning 
turns out to be incorrect. We consider bandlimited BMO functions in Section \5\ and establish some 
of their properties. In particular, we find that such a function / is either in L°° or that its samples 
{/(-)} are unbounded for any sampling rate s > 0. We exhibit a concrete example of such a func- 
tion, and study it in the context of our other results. 

We review some existing theory on bandlimited functions and the space BMO in Section [2j 
and discuss some preliminary results in Section [3] The main results of the paper are presented 
in Sections [4] and [5] We also develop our results for a class of general, nonuniformly spaced 
interpolation points, given by zeros of sine-type functions. The above discussion for uniformly 
spaced points is a special case. 

2 Background Material 

We will write f\ < f 2 if the inequality f\ < Cf 2 holds for a constant C independent of f\ and f 2 . We 
define f x > f 2 similarly, and write /1 ~ f 2 if both /i < f 2 and /1 > f 2 . For a set of points Y = {y k } 
and an extra element y, we denote the collection {jyt} \J{y} by Y, with | |F||/p := (| | Y\ \p + \y\ p ) l ' p 
and ||y||/~ :=max(||y||joo , |;y|). These conventions will be used throughout the paper. 

We first review a basic, alternative formulation of PW? , 1 < p <°°. An entire function / is said to 
be of exponential type b if 

£ = inf(/3:|/(z)|<e^UGC). 

We denote this by writing type(/) = b, and by type(/) = °° if b = 00 or / is not entire. By the 
Paley -Wiener-Schwartz theorem [9], PW£ can be equivalently described as the space of all entire 
functions with type(/) < %b whose restrictions to M. are in L p . It also follows that PW£ C PWS for 
p < q. Functions / G PW? satisfy the classical estimates H/'H^p < ^||/|Ilp and ||/(* + ic)|| i p < 
e %b\c\ ||y|| LP; respectively known as the Bernstein and Plancherel-Polya inequalities [fT0l[T2l . 

There is a rich and well-developed theory of nonuniform sampling for functions in PW£. We 
only cover a few aspects of it that we will need in this paper, and refer to lfT2ll and lfT5ll for more 
details. We consider a sequence of points X = {x^} C K, indexed so that x^ < x^+i. The separation 
constant of X is defined by A (X) = hr4 \xk+i —Xk\, and X is said to be separated if A (X) > 0. The 
generating Junction of X is given by the product 

Sx i: m TT / 1 *• 



sw^aa n (1--). <2) 

0<\x k \<r 



where 8x = 1 if G X and 8x = otherwise. For real and separated X, such a function S is said to 
be sine-type if the following conditions hold: 

(I) The product ^ converges and type(S) = nb < °°. 

(II) For any £ > 0, there are positive constants C\(e) and C2(e) such that whenever dist(z,X) > £, 

Ci(£)<e-^ Im ( z )l|5(z)|<C 2 (£). (3) 

It can be shown that condition (II) is equivalent to requiring that the bounds p} only hold in some 
half plane {z : |Im(z)| >c},c> 0. Furthermore, a sine-type function S also satisfies the bounds 
I £'(**) I ~ 1 an ^ f° rces X to satisfy sup^ \xk+\ — Xk\ < °° [10J. 



Now suppose the sequence X = {x{\ has a sine-type generating function S with type (5) = Tib. Let 
1 < p < oo. Then any / G ^W^ can be expressed in terms of its samples a^ = f(xk), 

/W = E g * «7 v Z) V (4) 



with uniform convergence on compact subsets of C. Conversely, for any {%} G l p , the series Q 



converges uniformly on compact subsets of C and defines a function / G PWF witha^ = /(x^j [10J. 



The simplest example of a sequence X with a sine-type generating function is the uniform 
sequence x^ = |, for which 5(z) = sm ^, ^ and the expansion (4) reduces to the WSK sampling 
theorem. More generally, any finite union of uniform sequences has a sine-type generating func- 
tion. As a more interesting example, the Bessel function 7o has real, separated zeros, satisfies 

Jo(z) =Jq(—z), and has the asymptotic formula Jq(z) = W^cos(z— f )(1 + 0(^)) as |z| — > °° and 

| argz| < n (see El). This implies that for sufficiently small £ > 0, S(z) = z/ (f )M^T^-) is a 
sine-type function with type (5) = %. Sequences X with sine-type generating functions are not the 
most general class for which / has an expansion of the form Q, but they have several convenient 
properties and cover some important cases encountered in applications, such as that of periodic 
interpolation points. Such sequences X and various properties of the series ([4]) have recently been 
studied in [4] in a computational context. 

The above results do not directly carry over to bounded functions / G PW^, but in this case we 
still have the following theorem j2|. 

Theorem. (Beurling) For a sequence X = {x^}, let N(X,I) be the number of Xk in an interval I. 
Then\\f\\ L „~\\f(X)\\ r forallfePWfifandonlyif 

, x • n N(X,\a,a + r)) 
D~(X) := limsup inf v ' L ^ > b. 

r— >oo a r 

D~(X) is called the lower uniform density of X. For a uniform sequence x^ = |, D~(X) = s, and 
Beurling's theorem implies that / G PW£° is uniquely determined by its samples if we oversample 
it beyond its Nyquist rate. 



We finally review a few properties of the Banach space BMO of functions with bounded mean 
oscillation, which has been studied extensively in connection with singular integral operators. It is 
defined by 

[/ = MWbmo = s f Ui I /(0 - j^T [f( s ) ds dt <oo\ 

where the supremumruns over all real intervals /. The quantity ||/||g M0 is technically a seminorm, 
since ||/||gMO = ll/ + c llsMO f° r an y constant c. Now for any g EL 1 , we denote its Hilbert trans- 
form by Jifg(z) := /"«, Jfl \ dt and its Riesz projections by 3? ± g := (g± iJifg) /2. We can then 
consider the "real" Hardy space H l (R), given by 



{/: 11/11* 



ri(») = WJ lli 1 + II^/IIli <0 ° 
Finally, it will also be useful to define the subspaces 

/CO 
f(t)dt = 0} 
-oo 

U 2 = {feH l (R):(\+t 2 )\&>+f(t)\eL°°} 

which are both norm dense in H l (R) [El[T3l. These spaces are all closely related, as the following 
theorem shows. 

Theorem. (Fefferman) BMO is the dual space ofH 1 (R). More specifically, we have the inequality 

1 



BMO ~ SU P 



geV \\8Wh 1 



f(t)g(t)dt 



where U can be taken as Ui or U2. Conversely, for any bounded linear functional L on H (R), 
?/zere 15 an f E BMO with \\L\\ ~ II/IIsmo- 

We write w = w + iv for the complex variable w in what follows. Let C ± = {w : ±v > 0} be the 
upper and lower half planes, and let P(w,t) = ^ , _>y2 + ,2 be the Poisson kernel on C + . Now define 
the square Q ar = {w : a < u < a + r,0 < v < r} . A measure /i on C + is said to be a Carleson 
measure if we have jY (\i) := sup f r "' ,a G R,r > 0] < °°. In other words, the measure /i 
of any square protruding from the real axis must be comparable to the length of its edge. The 
following theorem characterizes BMO in terms of such measures. 

Theorem. (Fefferman- Stein) Suppose f^ -jrvrdt < °° so that P{w, ■)■*•/ is well-defined. Then 



BMO 



t 2 +\' 
•* (v\V u , v (P(w,-)*f)\ 2 dudv^ ''" 



(5) 



A detailed discussion of BMO and the significance of these theorems can be found in [[81 or [fT3ll . 



3 Bandlimited Interpolation of Bounded Data 

In this section, we establish a preliminary result showing how adding an extra sample allows us to 
treat the bandlimited interpolation of bounded data, such as the noise model discussed in Section 
[TJ We define 

PW+ = { /entire : limsup / z~ 2 e- Kb \ Im ^f(z) \dz\ <°°\. (6) 

I ?--><*> J\z\=r ) 

The Plancherel-Polya inequality shows that PWu C PW&. Functions in PW& can be expanded in 
the following way. 

Theorem 1. Suppose X = {x k } C R is separated and has a sine-type generating function S with 
type(S) = nb, and let x^X. Iff e PW^andA = f(X), then 

_S(z) f, ,. S(zo) ( 1 1 \ 

/ z = a-^ + £ a* hm -^- , (7) 

S{x) k ~ x z ^zS'(x k ) \zo-x k x-x k J 

with uniform convergence of compact subsets of C Conversely, for any A G l°°, the series ([7|) 
converges uniformly on compact subsets ofC and f G PW^~. 

Proof. We use a standard complex variable argument. Assume z is in a closed ball B with z $■ X, 
and choose a real sequence {r n } with r n —> °° and dist({r n },X) > 0. We can then consider the 
integral 

2KlJ\ w \=r„ S(w) \z-w x-wj 
For sufficiently large n, it can be seen by calculating residues that 

J(r n 



\x k \<r n " l " u; V"' Xk X Xk 



\f(w)\e- 7tb ^ m ^ , 



r , . _S(z) _ S(z) / 1 

The inequalities Q and ([6]) imply that as r„ — V °° 

|J(r n )|<max|S(z)(z-x)|/* 

zes J\w\=r„ \w\ 

By letting z — > x k for each x k G 5, we obtain the formula ([7]) for all z G 5. For the other direction of 
Theoremfll we note that S has simple zeros at exactly X, so for z G K, |5(z) | < 2| |5'| |/,~dist(z,X). 
The Bernstein and Plancherel-Polya inequalities then show that for z G C and d = sup^ \x k +i —x k \ < 
°°. 

\S(z)\ < ||5|| L -min(dist(z,X), d)e Kb ^ m ^ . 

Now define the sets: 

/f = (LRe(w)J-min(l/2,A(X)),LRe(w)J+min(l/2,A(X))) 

7 2 = (-oo,L(Re( z )+i)/2j)\(/f|J/f) 
/ 3 = (L(Re(z)+i)/2j + l,oo)\(/j|J/f) 



Using the separation of X along with basic properties of lower Riemann sums, we have 
\f(z)\ < |aS(z)|+^ Im «l||A|| r £ min (dist(z,Z)^)|z-i| 



k= 



\Z-Xk\\X-Xk\ 



< 



7tb\im(z)\ | t , y \Xk+i - Xk\\z - x\ y \x k -x k -i\\z-x\ 



< ||A|Le» fc l Im «l(l+/ J^M 

< ||A|| r ^l Im «l(l+max(log|z|,0)), (8) 
which implies that / G PW^ + . □ 

This expansion can be compared with the series Q. It is essentially a nonuniform version of the 
classical Valiron interpolation formula considered in [3], in which the derivative of / at a point is 
used instead of the extra sample a, but the form considered here will be more convenient for our 
purposes. We also mention that the extra point x plays no special role in the collection X, and we 
isolate it mainly for notational convenience. If we pick any point Xj G X and let y^ = xi for k ^ j, 
yj = x and y = Xj, then Y = {y^} \Jy satisfies the conditions of Theorem IT] too. 



For any A G l°°, we call the function / given by (J7J) the bandlimited interpolation of A at X. Note 
that for any given 5,2 and X2 $ X, if g is the bandlimited interpolation of A U{^2} at X U{*2}> then 
g(z) — f(z) + cS(z) for some constant c. Moreover, if A G I 2 , then for any given x G" X we can 
always choose a so that / coincides with the series (|4]), or in the special case of uniformly spaced 
points X = {|}, the usual bandlimited interpolation given by the WSK series hb. 



We discuss an example of a PW^~ function that illustrates many of the typical properties of the 

m 



series (7). We use the uniform samples X = {k} and denote \j/(z) = Yf4, where T is the usual 



gamma function. The properties of \\r are discussed in depth in lfT4ll . 

Example. The function G\ (z) = sin(^z) y(— z) is in PW^\PW™ and satisfies a& = for k < and 
a k = (-!)% for k > 0. 



The function \y satisfies the estimate 



l|m fW = !, (9) 

~,|argz|<7r lOgZ 



so G\ is not bounded. With A = {%} given as above, Theorem [T] shows that for any x and a, the 
(unique) bandlimited interpolation of A at X is of the form G\(z) + csin(7Tz). It follows that the 
samples A have no bandlimited interpolation in PW™ . 

It will be instructive to isolate one property of G\ here. A classical formula of Gauss ([14], p. 240) 
shows that for integer k > 0, 

1 1 A" 1 1 2k ~ l 2 \ 

1 (k- 2 ) = G 1 (-k- 2 ) = (-lf(Z m+ i: m+ c\, (10) 

\m=l m=k I 



so as z —7- oo, |Gi(z)| grows logarithmically in between the integer samples. The same applies as 
Z — > — °°, even though the samples at k < are all zero. This can be interpreted as a nonlocal effect, 
where the sustained growth of |Gi| on the positive real axis, caused by the "bad behavior" of the 
samples at k > 0, induces growth on the negative real axis too. This property can be seen in the 
graph of G\ in Figure [1] It is also present in the bandlimited interpolation of Boche and Monich's 
example % = 0, k < 1, and % = (— \) k /\og(k+ 1), k > 1, where we take x = \ and a = 0. 





Figure 1: Left: The function G\(z). Right: The bandlimited interpolation of Boche and Monich's 
sequence. 

4 Bandlimited Interpolation of Random Data 

We can now state the main result of this paper. 

Theorem 2. Suppose X c M is separated and has a sine-type generating function S with type (5) < 
nb, and let x $■ X. Suppose also that A — {ak}\Ja is a collection of i.i.d. random variables 
uniformly distributed in [— a, a]. Let f be the bandlimited interpolation of A at X. Then almost 
surely, 

SUP^ f\f(t)\dt<oo. (11) 

r>0 lr J-r 

We make a few comments before proving Theorem[2] This result deals with the same situation 
discussed in Section [TJ even though it has been formulated slightly differently. In the notation of 
Section[T] we can take 7^ to be zero by linearity and only consider the noise N^. As we saw in Sec- 
tion [3j the extra sample a can be taken as deterministic and changed arbitrarily without affecting 
the result of Theorem [2j The exact probability distribution of A is also of little significance here, 
and the result holds more generally for any symmetric, finitely supported distribution. 

We split the proof of Theorem [2] into three lemmas for clarity. Our approach is to write the function 
/ as the sum of two parts, each with only zero samples in one direction along the real axis, and 
show that each one is almost surely bounded on that side. This shows directly that the nonlocal 
effect discussed in Section [3] does not occur. We then move to the deterministic setting and show 
that this one-sided boundedness forces a certain regularity upon the other side, resulting in the 



function having a bounded global average. 



For the rest of this section, we assume that X and S are as given in Theorem [2} without repeating 
the conditions on them every time. 

Lemma 3. For k such that Xk > 0, let {a^} be a collection of i.i.d. random variables uniformly 
distributed in [—a, a], let a^ = for all other k and let a = 0. Suppose f is the bandlimited 
interpolation of A atX. Then sup f<0 \f(t)\ < °° almost surely. 

Proof. We can assume that xq = min^ : Xk > 0) and x > 0, as the general case follows from the 



remarks after Theorem 



Let&yt 



a\ 



S'(xk)(x-Xk 



. Then we have 



£=0 







and the separation property shows that for some constant d, 

oo 2 oo 1 

var(^) - — £^77^2 



£var(&jt) = ^-E"E77 ,'/- 

*- 3 f n s'(x k y(x-x k y 



k=0' 

a 2 °° 

< —Y 

J k=0 



(dist(x,X) + X(X)\k-d\) 2 



By Kolmogorov's three-series theorem, Y,T=o^k converges almost surely. Now let 

f(t) £ a k ( 1 1 



(0 



Sit) 






o s ' ( x k) \t-x k x-x k 



It is easy to check that if £^L ^ converges, then lim f ^_ 00 ^(f) = YX=^k- Since |g(0) | < °°, it fol- 
lows by continuity that sup f<0 \g(t) | < °° almost surely. We also have sup f<0 \f(t) | < sup /<0 \g(t)\, 
which proves the lemma. □ 

Lemma 4. For any A £ l°°, let f be the bandlimited interpolation of A at X. Then for each c > 0, 



S(- + ic) 



< 



BMO 



Proof Applying Fefferman's duality theorem to the series ([7]) gives 



f(- + ic) 



S(-+ic) 



£ sup 

BMO heUi 



\hH 



Eakh(z) 



'k=- 



-'^S'ixk) \z + ic-x k x-x k 



dz 



Since h is finitely supported and the series ([7]) converges uniformly on compact sets, we can in- 
terchange the order of summation and integration. £P + h and £P~h are in L , so by analyticity we 



have 



H- + ic) 



S(- + ic) 



< 



BMO 



SU P Tuji 

heVi ll"llffi(R) 



sup 
heVi 



2%i 



y ak r ( ^ +h (z) + &~Kz) Hz) \ d 

.J^^S'ixk) J-oo \ z + ic-x k x-XkJ 

^i a.k&~h(xk — ic 



\h\u,\ 



< 



H l (R)k~- 
1 



z - sup 



\h 



h&i X ll"ll//l(R) fc=-oo 



S'{x k ) 



£ | £»"&(** -fc)|. 



Since X is separated, an elementary property of Hardy spaces ( ifTOll . p. 138) is that 

oo 

£ \&-h(x k -ic)\ < \\&-h\\ Li < \\h\\ HHR) , 

k=—oa 

which completes the proof. 



□ 



Lemma 5. For any A G l°°, let f be the bandlimited interpolation of A at X. Suppose that 
su Pt<o 1/(0 1 < °° and for some c> 0, f Jj'XZ\ e BMO. Then sup r>0 j- r f r _ r \f(t) \dt < oo. 



S{-+ic) 



Az+i) 



,Mi 



Proof. We assume c = 1 without loss of generality. Let f (z) = f(z)e z , g(z) — V( _ , 

sup r<0 |/(0I an d M-2 = sup f<0 |g(0l- The estimate (8) implies that f^-^rrr < °°, so |/ + | has a 
harmonic majorant on the upper half plane (see JH) and the reproducing formula f + (z) = P(z, •) * 
/ + holds for Im(z) > 0. We can then estimate 



sup|/ + (f + i)| < sup Mi / P(t + i,s)ds+ \f(s)\P(t + i,s)ds 
?<0 t<0 V J-™ Jo 



2 TtJo S 2 +l 



ds . 



This shows that M 2 < °°. Now for any fixed r > 0, 

;{t)\dt 
Y r ([\s{t)\dt - j_\g{t)\dt\ +M 2 



1 r 1 '•'■ 

-J \f(t+i)\dt < - 



2rJ- 



< 



< 



< 



2r 
1 

2~r 
1 '•'• 



rr rr rO rr 

/ \g(t)\dt- / \g(s)\ds + / |g(0l^- / \g(s)\ds 

JO J-r J-r J-r 



+ M 2 



K0 



2r 



< 2\\8\\BMO+ M 2- 



g(s)ds 



dt+M 2 



We finally use a Poisson integral again to move back to the real line. For Im(z) < 1, we have 
f-(z) = P(z-i, •)*/-(■+>')■ This gives 



2rJ- r uyn - 2rJ- r J-oc(t-s) 2 + l 



10 



1 f°° 

= e - — / |/(s + /)| (arctan (r + 5) + arctan(r — s)) ds 

\2rJ-2r JR\[-2r,2r] S l + 1 



Taking the estimate ([8]) into account again, we conclude that 



sup— / |/(0 \dt< 00. 
r>0 2r J-r 



We can now combine these lemmas to complete the proof. 



□ 



Proof of Theorem^ For any A e l°°, we can write the bandlimited interpolation / of A at X as 
f(z) = Mz)+f 2 (z) + |^, where h{x k ) = for ** < and f 2 {x k ) = for x k > 0. Applying 
Lemmas [3](5] on f\ (z) and fi{—z) and noting that S E L°° finishes the proof. □ 



The statistical incoherence in the samples A in Theorem^ is the reason we have the bounded 



average property (11), and it does not generally hold for bounded samples A. As an illustration of 
this, we return to the example function G\ from Section [3] and show that the average of |Gi(0l is 
unbounded. It suffices to consider t < 0. Let T be the tent function 




T{t) = { 2-2t ^<t< 1 . (12) 

otherwise 

It is clear that |sin(^)l > Ln=-°° T ( f + n ) > and the formula (J9j) implies that \yr(t)\ > ^ log |f| for 
sufficiently large t . This shows that 

n=1 

It follows that as r -> °°, \ f° r \ G\ (t) \dt>\ogr^°°. 

Figure [2] below shows an example of the bandlimited interpolation of random data. In the notation 
of Theorem 2l we use a realization of A with a = i, and take x k = k and x= \. We denote the 



resulting function by G2. The graphs in Figure [2] can be compared with the functions shown in 
Figure [T] in Section [3} Unlike those functions, it can be seen that G2 does not steadily grow over 
long time intervals. Intuitively, this shows how the effect of noisy samples on the bandlimited 
interpolation is in a sense well-controlled. 
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Figure 2: Left: The function G 2 (z) on [-100, 100]. Right: G 2 (z) on [-5000,5000]. 



5 Bandlimited BMO Functions 



In this section, we study some properties of bandlimited functions in the space BMO. Such func- 
tions have a somewhat different character than the examples we have seen so far. We fix a point c 
and define the space PW£ to be the following. 

pwt = if ■ typ e (/) < **> II/IIbm 0) c : = \f(c)\ + W/Wbmo < °°} 

The term |/(c) | resolves the ambiguity in the BMO seminorm for constant functions, and ||/|| BA /o c 
is a (full) norm. It will be shown below that the precise value of c is unimportant and that changing 
it gives an equivalent norm. Since / E BMO always satisfies /.Too'sxpfr < °° 181. the Paley- 
Wiener-Schwartz theorem implies that PW£ C PWu . We first give a version of the Plancherel- 
Polya inequality for PW£. 

Lemma 6. IffePWf, then ||/(- + *c)|| SM0 < ||/IU ^K 
Proof. The proof is similar to the PW£ case described in |fT2| . Define 



Rf(z) 



^(7tb+e)lm(z) 



TrS f{ - z+t) -U f(z+s)ds 



dt, 



< 



for complex z and real r. For each e > 0, R^ is a subharmonic function satisfying \Rj (z 

for z G M. and max(log \R^(z) |,0) — > as z — > i°°. Applying the Phragmen-Lindelof principle over 



BMO 



C+gives \R+(z + ic)\< 



BMO 



e (Kb+£)\c\ f or c > o, and we can repeat the argument with R e and 



C~ for c < 0. Taking the supremum over real z and r and letting e — > gives the inequality. □ 

We will now establish several basic properties of PW£. 

Theorem 7. Let f G PW?. Then the following statements hold. 

I: For each c Gl, /(• + *c) w uniformly Lipschitz continuous on K. 

Z7: For any fixed numbers c and c' , \f\sMO c ~ 11/llsMO c'- 

ZZ7: For any grven z G C, the point evaluation functional z — > f(z) is bounded on PW b *. 

IV: \\f\\ L . < \\f\\ BM0 . 
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Proof. We set b = 1 without loss of generality. We can prove all of the above statements by using 
the reproducing kernel-like function 



K(cJ) 



\c\ . (2tzN 
sin 1 



Ttt(t — c) \ c 



where c G M\{0} and N is any integer greater than \c\. As a function of t, K(c,t) is entire and 
satisfies 2% < type(K) < °°. For any / G PWf, 



f(t)K(c,t)dt=f(c)-f(0). 



This can be seen by observing that for T] = ±1, the function 7 ,, c _ c \ exp {2nir\Nz/c) f(z) has poles 
at c and with respective residues r\f{c) and — rj/(0). The estimation argument is very similar to 
the proof of Theorem [T] and we omit the details. 

We now suppose that / G PWf. We want to approximate the H 1 (R) norm of K(c, t) — K(c' ', ?), 
where c > 1 and c ; > 1. We first integrate the function % , — ^ -, c _, exp {2nir\Nz/c), where s G 
IR\{0, c}, and perform the same kind of calculation as before to find that 



Jf?K(c,s) 



1 



1 cexp(2niNs/c) cexp(—27CiNs/c) 



ns n(c — s) 2ns(c — s) 
c(cos(2nNs/c)- 1) 
7ts(c — s) 



2ns(c — s) 



Let A^ = max( \c] , \c'] ) and define the interval I w := [w — k,w + i]. We first consider the case 



where 1 < c < | and |c — c'| > i. Recalling that T is the tent function ( 12 ), we have 

\\K(c,-)-K(c',-)\\ Ll 

(2cT(2Nt/c + n) 2c'T(2Nt/c' + n) 



< 



r r, (2cT{2Nt 



16N 

< + 



< 



n 



\c — c \. 



2c 



+ ■ 



rdt + 



n\t{t-c') 



dt 



\(I°{JI c )7t\t(t-c)\ ' jR\{l {JI< J )7t\t(t-c' 



2d 



rdt 



Now suppose that 1 < c < | and \c — c'\ < k, so that N = 2. Some elementary estimates show that 



< 



< 



K(c,-)-K(c',-)\\ Li 

4 4 



1/2 
1/2 



c c 



5/2 

dt+ I max 

1/2 



4 sin(47rc/c') 



c n(c-c') 



4 sin(4^c'/c) 



c' n(c' — c) 



dt + 



\c — c \\t\ 



~ 3 / 2 dt 



\c — c 



\(-l/2,5/2) 

/i 
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Following the same arguments, we can also obtain the bound \\J$?K(c, •) — J4?K(d, •) || L i < |c — c'| 
for the above choices of c and d. By Fefferman's duality theorem and the fact that K(c, •) — 

K(d, •) G U2, we have 



BMO ~ 



1 



\K(c,-)-K(c',-)\\ H] 



(R) 



f(t)(K(c,t)-K(d,t))dt 



> 



f(c)-f(c f ) 



c — c 



(13) 



where the constant in the inequality is independent of c and d . Since the BMO seminorm is 
translation-invariant, the inequality (13) actually holds for all c,d E E. Combining this with 



Lemma k^ proves (I) and letting d — >■ c gives (IV). If we fix i? = |c — c'\ > 0, this also shows 



that ||/||bji/o + |/(c) I > |/(c')|. where the implied constant depends on 7?, and we can interchange 
c and d to get (II). Finally, the statement (III) is just (II) phrased in a different way. □ 

Remark. The closure of the set of uniformly continuous BMO functions under the BMO seminorm 
is called VMO, for vanishing mean oscillation. Theorem [7] (I) shows that PW£ C VMO. Note that 
there are two non-equivalent definitions of VMO in the literature, and we use the one given in [ 8 1 . 



Remark. Theorem [7] (IV) is a sharper form of the p = °° case of Bernstein's inequality. We mention 
that the opposite inequality does not generally hold (even if ||/|| SM0 is replaced by ||/||bmo c)> an( * 
there are functions / such that /' G PW^ but / £ PW£. 

Corollary 8. Let f G PW£. Then either f G PWf or there is no separated sequence X with 
D (X) > such that f(X) G f°. 

Proof. Suppose we have a separated X = {jCyt} with D~(X) > and f(X) G /°°. This means that 
for some large fixed r, every real interval / of length r contains a point x n EX. Theorem |7] (IV) 
then shows that for any t El, 



1/(01 



fix n ) + / f'(u)du 



\f(*n)\ + r 



BMO- 



□ 



Intuitively, Corollary ^1 says that an unbounded PW£ function is large in most places on the 
real line. It also shows that the bandlimited interpolation of bounded data A E l°° can never be 
in PWj* unless it is actually in PW^ . This occurs in spite of Lemma [4] and highlights a basic 
difference between PW£ and PW£, 1 < p < °°. In LemmaWj we generally cannot remove the fac- 



tor 



S(-+ic) 



from the inequality and conclude that / G BMO. In contrast, for A G l p , the series (4) 



can be used to find that 



f(-+'c) 
S{-+ic) 



E V (see 0H), which clearly implies /(• + ic) E IP and thus f EL P . 



We finally study an example of an unbounded PW£ function that illustrates the "largeness" 
property described above. 



Example. The function G 3 (z) = Y%= (-l) k s'm (fjfe\ is in PW* /3 \PW[ 



IV 
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To see this, we use the identity sinz = h [e lz — e~ lz ) to write G3 = G^ + + G3 , where P(w, •) * 
G3± = G-s±{w) for w G C ± , and then apply the Fefferman-Stein theorem ^ to each part. Let 
w = u + iv. We first note that by analyticity, 

\V(P(w r )*G 3 +)\ 2 = \VG 3+ (u + iv)\ 2 = 2\G' 3+ (w)\ 2 . 
for large k, the series defining G3 converges uniformly on compact sets, so 



Since 



sin 



nz 
3-2* 



we have 



< 



nz 



3-2* 



jV i 2v|G 3+ (w)| dudv 



,JV 2v 



d 1 



dwli 






42~ k w 



k=0 



dudv 



iit. 



< ,JV 2ve~~ 



< 2. 



K 



a=o 



3-2 k 



dudv 



Doing the same calculation with G3-, we find that G3 G ^^,3. On the other hand, G3 satisfies the 
identity G3(2z) = sin( 2 y I ) — G3(z). This implies that for integer n > 2, 

n- 1 / -jk \ 

G 3 (2») = (-l)>3(l) + I(-l)'^sin — 

k=o \ ^ / 



= (-1)" a(i; 



V3 



(n-2) 



so G3 ^ PWyy By Corollary 8 t the samples Gt,(X) are unbounded for any separated sequence X 
with D (X) > 0. It is interesting to note that such a function can still be bounded on a sequence X 
that is "very sparse" in the sense that D~(X) = 0. It is easy to check that G3P • 2 n ) = (-l) n G3(3) 
and G3(— z) = —Gj,{z), so G$(X) G /°° for the sequence x n — 3 • 2"sign(n). Some graphs of G3 are 
shown in Figure [3] below. 




5000 



Figure 3 : Left: The function G 3 (z) on [- 100, 100] . Right: The absolute value of G 3 (z) on [0, 5000] . 
The peaks at powers of 2 are clearly visible, as well as a self-similarity effect at different scales. 
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